PTEP 



Prog. Theor. Exp. Phys. 2012, 00000 (25 pages) 
DOI: 10.1093/ptep/OOOOOOOOOO 



Construction of gauge-invariant variables of 
linear metric perturbations on an arbitrary 
background spacetime 

Kouji Nakamura^ 

^ TAMA Project, Optical and Infrared Astronomy Division, 
National Astronomical Observatory of Japan, 
2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 
' * E-mail: kouji.nakamura@nao.ac.jp 

(N ■ 



An outline of a proof of the decomposition of linear metric perturbations into gauge- 
invariant and gauge-variant parts on the an arbitrary background spacetime which 
admits ADM decomposition is discussed. We explicitly construct the gauge-invariant 
and gauge-variant parts of the linear metric perturbations through the assumption of 
the existence of some Green functions. We also confirm the result through another 
approach. This implies that we can develop the higher-order gauge-invariant perturba- 
tion theory on an arbitrary background spacetime. Remaining issues to complete the 
^ , general-framework of the general-relativistic higher-order gauge-invariant perturbation 

bX)' theories are also discussed. 
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2 ■ !• Introduction 

'NT 

ly-^ I Perturbation theories are powerful techniques in many area of physics and the develop- 

O ■ ments of perturbation theories lead physically fruitful results and interpretations of natural 

phenomena. 

In physics, researchers want to describe realistic situations in a compact manner. Exact 
. solutions in a theory for natural phenomena are candidates which can describe realistic 

^ ! situations. However, in many cases, realistic situations are too complicated and often difficult 

to describe by an exact solution of a theory. This difficulty may be due to the fact that exact 
solutions only describe special cases even if the theory is appropriate to describe the natural 
phenomena, or may be due to the lack of the applicability of the theory itself. Even in the 
case where an exact solution of a theory well describes a physical situation, the properties 
of the natural system will not be completely described only through the exact solution. In 
natural phenomena, there always exist "fluctuations" . In this case, perturbative treatments 
of the theory is a powerful tool and researchers investigate perturbative approach within a 
theory to clarify the properties of fluctuations. 

General relativity is one of theories in which the construction of exact solutions is not so 
easy. Although there are many exact solutions to the Einstein equation [1], these are often 
too idealized. Of course, there are some exact solutions to the Einstein equation which well- 
describe our universe, or gravitational field of stars and black holes. These exact solutions by 
itself do not describe fluctuations around these exact solutions. To describe them, we have 
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to consider the perturbations. Therefore, general relativistic linear perturbation theory is a 
useful technique to investigate the properties of fluctuations around exact solutions [2-4]. 
Through these linear perturbation theories, we can describe fluctuations such as the density 
or the temperature fluctuations of our universe, gravitational waves from self-gravitating 
objects. 

On the other hand, higher-order general-relativistic perturbations also have very wide 
applications. Among these applications, second-order cosmological perturbations are topical 
subject [5-9] due to the precise measurements in recent cosmology [10]. Higher-order black 
hole perturbations are also discussed in some literature [11]. Moreover, as a special example 
of higher-order perturbation theory, there are researches on perturbations of a spherical 
star [12], which arc motivated by the investigation of the oscillatory behaviors of a rotating 
neutron star. Thus, there are many physical situations to which general relativistic higher- 
order perturbation theory should be applied. 

As well-known, general relativity is based on the concept of general covariance. Intuitively 
speaking, the principle of general covariance states that there is no preferred coordinate 
system in nature, and the notion of general covariance is mathematically included in the 
definition of a spacctime manifold in a trivial way. This is based on the philosophy that 
coordinate systems are originally chosen by us, and that natural phenomena have nothing 
to do with our coordinate system. Due to this general covariance, the "gauge degree of free- 
dom", which is unphysical degree of freedom of perturbations, arises in general-relativistic 
perturbations. To obtain physical results, we have to fix this gauge degrees of freedom or 
to extract some invariant quantities of perturbations. This situation becomes more com- 
plicated in higher-order perturbation theory. In some linear perturbation theories on some 
background spacetimes, there are so-called gauge-invariant perturbation theories. In these 
theories, one may treat only variables which are independent of gauge degree of freedom with- 
out any gauge fixing. Therefore, it is worthwhile to investigate higher-order gauge-invariant 
perturbation theory from a general point of view. 

According to these motivations, the general framework of higher-order general-relativistic 
gauge-invariant perturbation theory has been discussed in some papers [13, 14] by the 
present author. Although these development of higher-order perturbation theory was origi- 
nally motivated by researches on the oscillatory behaviors of a self-gravitating Nambu-Goto 
membrane [15], these works are applicable to cosmological perturbations and we clarified 
the gauge-invariance of the second-order perturbations of the Einstein equations [9, 16, 17]. 

In Ref. [13], we proposed a procedure to find gauge- invariant variables for higher-order 
perturbations on an arbitrary background spacetime. This proposal is based on the sin- 
gle assumption that we already know the procedure to find gauge-invariant variables for 
linear-order metric perturbations (Conjecture 1 in Sec. 2 in this paper). Under the same 
assumption, we summarize some formulae for the second-order perturbations of the curva- 
tures and energy-momentum tensor for the matter fields in Refs. [14, 16]. Confirming that 
the above assumption in the case of cosmological perturbations is correct, in Refs. [9], we 
develop the second-order gauge-invariant cosmological perturbation theory. Through these 
works, we find that our general framework of higher-order gauge-invariant perturbation the- 
ory is well-defined except for the above assumption for linear-order metric perturbations. 
Therefore, we proposed the above assumption as a conjecture in Ref. [16]. If this conjecture is 
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true, the higher-order general-relativistic gauge-invariant perturbation theory is completely 
formulated on an arbitrary background spacetime and has very wide applications. 

The main purpose of this paper is to give a scenario of a proof of this conjecture based 
on the premise that the background spacetime admits ADM decomposition. We explicitly 
construct the gauge-invariant and gauge-variant parts of the linear metric perturbation. 
Although some special modes are excluded in the proof in this paper, we may say that the 
above conjecture is almost correct for linear-order perturbations on an arbitrary background 
spacetime. This paper is the full paper version of our previous short letter [18]. 

The organization of this paper is as follows. In Sec. 2, we briefly review the general frame- 
work of the second-order gauge-invariant perturbation theory developed in Rcfs. [13, 14] 
and the above conjecture is also declared as Conjecture 1 in this review. In Sec. 3, we 
give a scenario of a proof of Conjecture 1. We note that we assume that the existence of 
Green functions for two elliptic differential operators in our outline of a proof. Therefore, 
the modes which belong to the kernel of these two elliptic differential operators are excluded 
in this paper. Since we use tricky logic in our outline in Sec. 3, we reconsider the deriva- 
tion of gauge-transformation rules through an alternative approaches in Sec. 4 to check the 
consistency of our result in Sec. 3. The final section (Sec. 5) is devoted to summary and dis- 
cussions. Throughout this paper, we use the covariant decomposition of symmetric tensors 
on Riemannian manifold which was developed by York [26]. We review York's discussions on 
this covariant decomposition in Appendix A. Although his discussions are for tensors on a 
closed manifold, we use his decomposition on a finite region in a Riemannian manifold with 
boundaries assuming the existence of Green functions for two elliptic differential operators. 

We employ the notation in Refs. [13, 14] and use abstract index notation [19]. We also 
employ natural units in which the velocity of light satisfies c = 1. 

2. General framework of the higher-order gauge-invariant perturbation theory 

In this section, we briefly review the general framework of the gauge-invariant perturbation 
theory developed in Ref. [13]. The aim of this section is to emphasize that Conjecture 1 is 
an important premise of our general framework. In Sec. 2.1, we review the notion of the 
gauge in general relativity and gauge degree of freedom in general-relativistic perturbations. 
In Sec. 2.2, the definition of perturbations in general relativity and its gauge-transformation 
rules are reviewed. When we consider perturbations in any theory with general covariance, 
we have to exclude these gauge degrees of freedom in the perturbations. To accomplish this, 
gauge-invariant variables of perturbations are useful, and these are regarded as physically 
meaningful quantities. In Sec. 2.3, a procedure to find gauge-invariant variables of pertur- 
bations is explained, which was developed in Ref. [13]. We emphasize that the ingredients of 
this section do not depend on the details of the background spacetime, except for Conjecture 
1. 

2.1. Notion of "gauge" in general relativity 

General relativity is a theory based on general covariance, which intuitively states that there 
is no preferred coordinate system in nature. Due to this general covariance, the notion of 
"gauge" is introduced in the theory. Sachs [20] is the first person who pointed out that there 
are two kinds of "gauges" in general relativity. He called these two "gauges" as the first- 
and the second-kind of gauges, respectively. The distinction of these two different notion 
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of "gauges" is an important premise in the arguments in Sec. 3. Therefore, first of all, we 
remind the difference of these two "gauges" . 

The first kind gauge is a coordinate system on a single manifold M- On a manifold, we can 
always introduce a coordinate system as a diffeomorphism ipa which maps from an open set 
Oa (Z M to tpaiOa) C M"+^ (whcrc n + 1 = dimA^). This coordinate system ipa is called 
gauge choice of the first kind. If we consider another open set C Al, we have another 
gauge choice -.O^^ ipp{0^) C M""'"-'^. If Oq n 7^ 0, we can consider the diffeomorphism 
V'/? ° i'a^-i which is a coordinate transformation : tpa{Oa H O^) C M"+^ i-)- '^^{0^ fl O^) C 
]^n+i fpj^jg coordinate transformation il)^ ° i/^a^ is also called gauge transformation of the 
first kind in general relativity. 

On the other hand, the second kind gauge appears in perturbation theories in a theory 
with general covariance. This is the main issue of this paper. In perturbation theories, 
we always treat two spacetime manifolds. One is the physical spacetime M which is our 
nature itself and we want to clarify the properties of M. through perturbations. Another 
is the background spacetime Mo which has nothing to do with our nature but is prepared 
by hand for perturbative analyses. Let us denote the physical spacetime by {M,Q) and 
the background spacetime by (MqjQq), where Q is the collection of tensor fields on M, 
and Qq is the collection of the background values on A4o for the collection Q on A^. The 
gauge choice of the second kind is the point identification map X : Mo M [20, 21]. We 
have to note that the correspondence X between points on Mo and M is not unique to 
the perturbation theory with general covariance. General covariance intuitively means that 
there is no preferred coordinate system in the theory. Due to general covariance, we have 
no guiding principle to choose the identification map X. Actually, as a gauge choice of the 
second kind, we may choose a different point identification map y from X. This implies 
that there is degree of freedom in the gauge choice of the second kind. This is the gauge 
degree of freedom of the second kind in perturbation theory. In this understanding, the gauge 
transformation of the second kind is a change X y oi the identification map A^o '-^ M- 

2.2. Perturbations in general relativity 

To formulate the above second-kind gauge in more detail, we introduce an infinitesimal 
parameter A for perturbations. Further, we consider the (n + 1) + 1-dimensional manifold 
A/" = Al X R, where n + 1 = dim A^ and A G R. The background spacetime Mq = A^I^^q and 
the physical spacetime M = Mx = AA|j^^^ are also submanifolds embedded in the extended 
manifold A". Each point on Af is identified by a pair, (p, A), where p € Mx, and each point 
in the background spacetime Mq in A" is identified by A = 0. Through this construction, 
the manifold J\f is foliated by (n + l)-dimensional submanifolds Mx of each A, and these 
are diffeomorphic to the physical spacetime M and the background spacetime Mo- The 
manifold A/" has a natural differentiable structure consisting of the direct product of M 
and M. Further, the perturbed spacetimes A^a for each A must have the same differential 
structure by this construction. 

If a tensor field Qx is given on each Mx, Q\ is automatically extended to a tensor field 
on J\f by Q{p, A) := Qx{p), where p G Mx- Tensor fields on J\f obtained through this con- 
struction are necessarily "tangent" to each Mx, i-e., their normal component to each Mx 
in A/" identically vanishes. To consider the basis of the tangent space of A/", we introduce 
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the normal form {dX)a and its dual {d/dX)"", which are normal to each A4\ in Af. These 
satisfy {dX)a (d/dX)"' = 1. {dX)a and {d/dXY are normal to any tensor field extended from 
the tangent space on each Mx through the above construction. The set consisting of (dA)a, 
{d/dXY, and the basis of the tangent space on each M.\ is regarded as the basis of the 
tangent space of N . 

To define perturbations of an arbitrary tensor field (J, we have to compare Q on the 
physical spacetime M\ with Qq on the background spacetime Mq through the introduction 
of a gauge choice of the second kind. The gauge choice of the second kind is made by 
assigning a diffeomorphism Xx ■ M ^ N such that Xx : Mq — > Mx- The pull-back X^, 
which is induced by the map Xx-, maps a tensor field Q on M.^ to a tensor field X^Q on Mq. 
Once the definition of the pull-back of the gauge choice Xx is given, the perturbations of a 
tensor field Q under the gauge choice Xx are simply defined by the evaluation of the Taylor 
expansion at A^o in A/": 

:= XlQx\j^^ = Qo + A^i^Q + ^A^^Jq + 0(A3), (1) 

where '^^^^Q and ^"^^Q are the first- and the second-order perturbations of Q, respectively. 

We also note that these perturbations completely depend on the gauge choice Xx- When 
we have two different gauge choices Xx and 3^;^, wc have two different representations of the 
perturbative expansion of the pulled-backed variables '^x^Amo ""^ '^A^-'^I^Uo' 

■■= 3^IQAU„ = Qo + A5Q + ^A2(Jg + 0(A3), (2) 

Although these two representations of the perturbations are different from each other, these 
should be equivalent because of general covariance. This equivalence is guaranteed by the 
gauge-transformation rules between two different gauge choices. The change of the gauge 
choice from Xx to 3^^ is represented by the diffeomorphism <tx '■= (-^a)"^ ° ^A- This diffeo- 
morphism $A is the map $a : M-o M-o for each value of A G R and does change the point 
identification. Therefore, the diffeomorphism ^x is regarded as the gauge transformation ^x 
: Xx yx- The gauge transformation ^x induces a pull-back from the representation "^Qx of 
the perturbed tensor field Q in the gauge choice Xx to the representation -^Qx in the gauge 
choice yx- Actually, the tensor fields '^Qx and -^Qx, which are defined on A^o, are connected 
by the linear map as ■^Qx = ^x'^Qx- According to generic arguments concerning the 
Taylor expansion of the pull-back of tensor fields on the same manifold [7, 8], it should be 
expressed the gauge transformation ^l^Qx in the form 

yq = c&l^Q = + A£^,/Q + y {i:^,,, + 4 J -^g + 0(A3), (3) 

where the vector fields S,^^-^ and ,^^2) are the generators of the gauge transformation ^x- 
Substituting Eqs. (1) and (2) into Eq. (3), we obtain the gauge-transformation rules for 
perturbations ^^Q and ^'^Q as follows: 

^'JQ-^'Jq = £^,,,Qo, (4) 

%^Q-^'Jq = 2£^J'Jq + {£^^,^+£1^]Qo- (5) 

The notion of gauge invariance considered in this paper is the order by order gauge 
invariance proposed in Ref. [16]. We call the fcth-order perturbation is gauge invariant 
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iff 

= ^y^Q (6) 

for any gauge choice Xx and yx. Through this concept of the order by order gauge invari- 
ance, we can decompose any perturbation of Q into the gauge-invariant and gauge-variant 
parts, as shown in Ref. [13]. In terms of these gauge-invariant variables, we can develop the 
gauge-invariant perturbation theory. However, this development is based on a non-trivial 
conjecture, i.e.. Conjecture 1 for linear order metric perturbations as explained below. 

2.3. Gauge-invariant variables 

Inspecting the gauge-transformation rules (4) and (5), we define gauge-invariant variables 
for metric perturbations and for perturbations of an arbitrary matter field. First, we consider 
the metric perturbation. The metric gab on M, which is pulled back to Mo using a gauge 
choice X\, is expanded in the form of Eq. (1): 

Xl9ah = 9ab + Xxhab + Y^<^b + 0^{X), (7) 

where gab is the metric on A4o- Of course, the expansion (7) of the metric depends entirely 
on the gauge choice X\. Nevertheless, henceforth, we do not explicitly express the index of 
the gauge choice Xx if there is no possibility of confusion. 

Based on these setups, in Ref. [13], we proposed a procedure to construct gauge- invariant 
variables for higher-order perturbations. Our starting point to construct gauge-invariant 
variables is the following conjecture for the linear-order metric perturbation hab defined by 
Eq. (7): 

Conjecture 1. // there is a symmetric tensor field hab of the second rank, whose gauge 
transformation rule is 

2^ab — xhab = £(,(i)gabi (8) 

then there exist a tensor field Hab O'^d a vector field X"' such that hab is decomposed as 

hab ='■ T-Lab + £xgab, (9) 

where Hab o.f-d X"" are transformed as 

yH-ab - xHab = 0, yX" - xX'^ = 4(1) (10) 

under the gauge transformation (8), respectively. 

In this conjecture, T-Lab is gauge-invariant and we call T-Lab a-s gauge-invariant part of the 
perturbation hab- On the other hand, the vector field X'^ in Eq. (9) is gauge dependent, and 
we call X" as gauge-variant part of the perturbation hab- 

The main purpose of this paper is to give an outline of a proof of Conjecture 1. In the case 
of cosmological perturbations on a homogeneous and isotropic universe, we confirmed that 
Conjecture 1 is almost correct except for some special modes of perturbations, and then 
we could develop the second-order cosmological perturbation theory in a gauge-invariant 
manner [9]. On the other hand, in the case of the perturbation theory on an arbitrary 
background spacetime, this conjecture is a highly non-trivial statement due to the non- 
trivial curvature of the background spacetime, though its inverse statement is trivial. We 
will see this situation in detail in Sec. 3. 
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Before going to our outline of a proof of Conjecture 1, we explain how the higher-order 
gauge-invariant perturbation theory is developed based on this conjecture, here. Through 
this explanation, we emphasize the importance of Conjecture 1. 

As shown in Ref. [13], the second-order metric perturbations lab are decomposed into 
gauge-invariant and gauge- variant parts through Conjecture 1. Actually, using the gauge- 
variant part X"' of the linear-order metric perturbation hab, we consider the tensor field Lab 
defined by 

Lab '■= Lb - '^£xhab + £x9ab- (H) 

Through the gauge-transformation rules (5) and (10) for lab £^nd X", respectively, the gauge- 
transformation rule for this variable Lab is given by 

yLab - xLab = £a9ab, := C(2) + [^(1)' • (1^) 

This is identical to the gauge- transformation rule (8) in Conjecture 1. Therefore, we may 
apply Conjecture 1 to the variable Lab and we can decompose it as 

Lab = ^ah + ^YOah, (13) 

where the gauge-transformation rules for Cab a^id are given by 

yCab - xJCab = 0, yF'^-;ti^« = ^^2) + [C(i),Xf. (14) 

Thus, we have accomplished the decomposition of the second-order metric perturbation lab 
as 

lab ='■ C,ab + '^£xhab + {£y — £'x) Sab- (15) 

Furthermore, as shown in Ref. [13], using the first- and second-order gauge- variant parts, 
X" and Y"-, of the metric perturbations, the gauge-invariant variables for an arbitrary tensor 
field Q other than the metric are given by 

(i)q := ^'^Q-£xQo, (16) 

:= (^)Q-2£x^'^Q-{£y-£'x}Qo. (17) 

It is straightforward to confirm that the variables ^^^Q and ^"^^Q defined by (16) and (17), 
respectively, are gauge invariant under the gauge-transformation rules (4) and (5), respec- 
tively. We have to emphasize that not only gauge-invariant parts of the metric perturbations 
but also gauge-variant parts X" and for metric perturbations play crucial role in these 
systematic construction of gauge invariant variables ^^^Q and ^^^Q through Eqs. (16) and 
(17). 

Equations (16) and (17) have an important implication. To see this, we represent these 
equations as 

Wq = (^^Q + £xQo, (18) 

= (■'^Q + 2£x^^^Q + {£y-£^Qo. (19) 

These equations imply that any perturbation of first and second order can always be decom- 
posed into their gauge-invariant and gauge-variant parts as Eqs. (18) and (19), respectively. 
The decomposition formulae (18) and (19) are important consequences in our higher-order 
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gauge- invariant perturbation theory. Actually, in Ref. [14], we have derived the formulae for 
the perturbations of the spacetime curvatures and showed that all of these are decomposed as 
Eqs. (18) and (19). In addition to the spacetime curvatures, in Ref. [16], we also summarized 
the formulae for the perturbations of the energy-momentum tensors for a single scalar field, 
a perfect fluid, and an imperfect fluid, and showed that all these energy- momentum tensors 
and the equations of motion for the matter fields are decomposed into their gauge-invariant 
and gauge- variant parts as Eqs. (18) and (19). As a result of these decompositions, we can 
easily show that order by order perturbative equations for any equation on an arbitrary back- 
ground spacetime are automatically given in gauge-invariant form [14, 16]. Furthermore, we 
explicitly derived of the second-order Einstein equations for cosmological perturbations in a 
gauge-invariant manner [9, 16]. 

We can also expect that the similar structure of equations of the systems will be main- 
tained in any order perturbations and our general framework will be applicable to any order 
general-relativistic perturbations. Actually, decomposition formulae for the third-order per- 
turbations in two-parameter case, which correspond to Eqs. (18) and (19), are explicitly given 
in Ref. [13]. Therefore, the similar development is possible for the third-order perturbations. 
Since we could not find any difficulties to extend higher-order perturbations [13] except 
for the necessity of long cumbersome calculations, we can construct any order perturbation 
theory in gauge-invariant manner, recursively. 

We note that, through the above recursive procedure, we can find gauge-invariant variables 
for any perturbative variables without any gauge fixing. The concept of gauge invariance of 
perturbations should be equivalent to "complete gauge-fixing" . Therefore, we may say that 
our procedure gives a complete gauge-fixing without any explicit gauge-fixing. We also note 
that the specification of gauge-variant parts is not unique. The different specifications of 
gauge-variant variables X" and correspond to the different gauge-fixing as explicitly 
shown in Ref. [24]. In many literature, the explicit gauge-fixing procedures were proposed 
and these correspond to the explicit specification of the gauge- variant parts X" and ¥"■. 
However, in this paper, we do not explicitly specify these gauge-variant variables, though we 
can specify these gauge-variant variables at any time. This is due to the fact that the key 
idea of our recursive procedure to find gauge-invariant variables is not in the explicit form 
of the gauge-variant parts but in the gauge-transformation rules of gauge-variant variables 
of metric perturbations. For example, if we explicitly specify the first-order gauge-variant 
part X"', it will be difficult to construct gauge-invariant variables for the second-order metric 
perturbations, because the gauge-transformation rule of the gauge- variant variable X" is used 
to find gauge-invariant variables for the second-order metric perturbations in our recursive 
procedure. 

Finally, we have to emphasize that the above general framework of the higher-order gauge- 
invariant perturbation theory are independent of the explicit form of the background metric 
gab except for Conjecture 1, and are valid not only in cosmological perturbation case but also 
the other generic situations if Conjecture 1 is true. This implies that if we prove Conjecture 
1 for an arbitrary background spacetime in some sense, the above general framework is 
applicable to perturbation theories on any background spacetime. This is the reason why 
we proposed Conjecture 1 in Ref. [16]. In the next section, we give a scenario of a proof of 
Conjecture 1 on an arbitrary background spacetime which admits ADM decomposition. 
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3. Decomposition of the linear-order metric perturbation 

Here, we give a scenario of a proof of Conjecture 1 in Sec. 3.1 through the assumption of the 
existence of some Green functions of elliptic type derivative operators. This scenario is just 
an extension of the proof in the case of cosmological perturbations. The comparison with 
the case of cosmological perturbations in Refs. [9] is discussed in Sec. 3.2. 

3.1. A scenario of a proof of Conjecture 1 

Now, we give a scenario of a proof of Conjecture 1 on an arbitrary background spacctime. To 
do this, we assume that the background spacetimes admit ADM decomposition. Therefore, 
the background spacetime AAq (at least the portion of A^o which we are addressing) con- 
sidered here is n -|- 1-dimensional spacetime which is described by the direct product R x S. 
Here, R is a time direction and S is the spacelike hypersurface (dim E = n) embedded in 
A^o- This means that M.q is foliated by the one-parameter family of spacelike hypersurface 
S(t), where t G M is a time function. In this setup, the metric on Mq is described by 

gab = -a\dt)a{dt)b + qij{dx' + p'dt)a{dx^ +I3^dt)b, (20) 

where a is the lapse function, is the shift vector, and qab = Qij{dx'^)a{dx'^)i, is the metric 
on S(t). 

Since the ADM decomposition (20) of the metric is a local decomposition, we may regard 
that the arguments in this paper are restricted to that for a single patch in Mo which is 
covered by the metric (20). Further, we may change the region which is covered by the 
metric (20) through the choice of the lapse function a and the shift vector /5*. The choice of 
a and /3* is regarded as the first-kind gauge choice explained in Sec. 2.1, which have nothing 
to do with the second-kind gauge as emphasized in Sec. 2.1. Since we may regard that the 
representation (20) of the background metric is that on a single patch in Mo, in general 
situation, each S may have its boundary 9S. For example, in asymptotically fiat spacetimes, 
dT, includes asymptotically flat regions [19]. Furthermore, if necessary, we may regard that 
S(t) is a portion of the spacelike hypersurface in Mo and add disjoint components to the 
boundary 9S. For example, when the formation of black holes occurs, we may exclude the 
region inside the black holes from S. In any case, when we consider the spacelike hypersurface 
S with boundary dT,, we have to impose appropriate boundary conditions at the boundary 
5E. 

To consider the decomposition (9) of hab, we first consider the components of the metric 

hab as 

hab = htt{dt)a{dt)b + 2hti{dt)(^a{dx% + hij{dx\{dx^)b. (21) 

The components hu, ha, and hij are regarded as a scalar function, components of a vector 
field, and the components of a symmetric tensor field on the spacelike hypersurface S, 
respectively. Under the gauge-transformation rule (8), the components {hu, hu, hij} are 
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transformed as 

yhtt-xhtt = 2dt^t--idta + (3'Dia-l3^f5'Kij)Ct 

a \ 

WD'a - aP^D'Pk - ^'^^Djo^ ^i, (22) 
yhti-xhti = dtCi + Di^t-l{Dia-p^Kij)^t-lM,^^j, (23) 

yhij-xhij = 2D^iCj) + lKij^t-^P^Kij^k, (24) 
where M- ■' is defined by 

mJ:= -a'^K^. + P^p'^Kki - pWia + aDiP^ . (25) 

Here, Kij is the components of the extrinsic curvature of S in Mo and Di is the covariant 
derivative associate with the metric qij {Diqji^ = 0). The extrinsic curvature Kij and its trace 
K are related to the time derivative of the metric qij by 

1 



Kii =- 



d 

—qij - DiPj - DjPi 



K := q'^Kij. (26) 



We also note that the gauge-transformation rules (22)-(24) represent the gauge- 
transformation of the second kind which have nothing to do with the gauge-degree of freedom 
of the first kind as explained in Sec. 2.1. We have to emphasize that the main purpose of this 
paper is to show how to exclude this gauge degree of freedom of the second kind inspecting 
gauge-transformation rules (22)-(24). 

The essence of our strategy for the explicit construction of the gauge-invariant and gauge- 
variant parts of the linear metric perturbation is already given in our short paper [18]. In 
Ref. [18], we consider the simple case where a = 1 and = 0. Our strategy is as follows: we 
first assume that the existence of the variables Xt and Xi whose gauge-transformation rules 
are given by 

yXt - xXt = 6, (27) 

yXi — xXi = ^i. (28) 

This assumption is confirmed through the explicit construction of the gauge- variant part of 
the linear-order metric perturbation below. Similar technique is given by Pereira et al. [22] 

in the perturbations on Bianchi type I cosmology. 

Inspecting gauge-transformation rules (22)-(24), we define the symmetric tensor field Hab 
whose components are given by 

Htt ■■= htt + -{dta + pWia-P^p'Kij)Xt 
a 

a \ 

WD^a - aP'^D^Pk - P'P^DjO^ Xi, (29) 

Hti := hti + -{Dia-l3^Kij)Xt + -M.^Xj, (30) 
2 2 

Hij := hij KijXt -\ — P KijXk. (31) 
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The gauge transformation rules (22)-(24) and our assumptions (27) and (28) give the gauge- 
transformation rules of the components of Hab as follows: 

yHtt-xHtt = 2dtiu (32) 
yHti - xHti = dtii + (33) 
yHij — xHij = 2£)(j^j) . (34) 

Since the components Ha and Hij are regarded as componets of a vector and a symmetric 
tensor on S(t), respectively, we may apply York's decomposition reviewed in Appendix A 
to Hti and Hif. 

Hti =: Dih^vL) + h(v)i, D^h(y)i = 0, (35) 
Hij =: ^qijh(^L) +h(T)ij, Q^^h^T)ij=0, (36) 

^{T)ij =• {Lh(jv))^j+h(TT)ij^ D''h(xT)ij (37) 

where {Lh(TV))ij is defined by [see Eq. (A5) in Appendix A] 

2 I 

{Lh(TV))ij '■= Dih{TV)3 + Djh(Tv)i - ~^i3^ h{TV)l- (38) 

Equations (33) and (34) give the gauge-transformation rules for the variables h(^vL): 

h{L), h(T)ij, h(TV)ij and h(TT)ij- 
First, we consider the gauge-transformation rule (33) in terms of the decomposition (35): 

yHti — xHti = Di — xh(yL)) + {:i^{v)i — xh(y)i) = dtii + (39) 

Taking the divergence of this gauge-transformation rule and through the property D^h(y)i = 
0, we obtain 

A i^ivL) - xhivL)) = D'dtii + A(t. (40) 

In this paper, we assume that the existence of the Green function A^^ of the Laplacian 
A := D^Di. Then, we easily obtain the gauge-transformation rule for h(yi,-^ as 

y\yL) - X\VL) = 6 + A-^D^dtik. (41) 

where we have ignored the modes which belong to the kernel of the derivative operator A. 
Substituting Eq. (41) into Eq. (39) we obtain the gauge-transformation rule for the variable 

yhiy)i-xh^v)i = dt^i-DiA-^D^'dtCk. (42) 

The gauge-transformation rules for /i^j^) and /i(T)ij are given from Eq. (34). Since we 
consider the decomposition (36), the gauge-transformation rule (34) is given by 

yHij - xHij = -Qij {:^{L) - xh(^L)) + (yi(r)y - xh(j)ij) = ^D{i^j)- (43) 
Taking the trace of Eq. (43), we obtain 

yhiD-xh^L) = 'iD^i. (44) 
The traceless part of Eq. (43) is given by 

yh{T)ij - x\T)ij = m)-- . (45) 
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Note that the variable /i(T)jj is also decomposed as Eq. (37) and the gauge-transformation 
rules for the variable ^(t),j is given by 

yh{T)ij - xh(^T)ij = {L (yi(TV) - xh{TV)))^j + yh{TT)ij - xh(TT)ij 

= {LOij- (46) 
Taking the divergence of Eq. (46) , we obtain 

'D'^ {yh(TV)i - xhiTV)i - 6) = 0, (47) 
where the derivative operator T>^^ is defined by 

qij^ + ^1 _ Dij^J + jiiJ_ (48) 

Here, R^^ is the Ricci curvature on S. Properties of the derivative operator 2?*-' are discussed 
in Appendix A. Here, we assume the existence of the Green function of the derivative operator 
"D*^ and ignore the modes which belong to the kernel of the derivative operator V^^ . Then, 
we obtain 

yh(TV)l - xh(Ty)i = ^l. (49) 

Substituting Eq. (49) into (46), we obtain 

yh{TT)ij - xh(TT)ij = 0. (50) 

In summary, we have obtained the gauge-transformation rules for the variables Hu, h(yL)j 
h{v)i, h(^L), h(Tv)i, and /t(rr)y as follows: 

yHtt-xHtt = 2dt^t, (51) 

yh^VL) - xh^VL) = Ct + A~^D%^k, (52) 

yhiv)^ - xh^v)i = dtii-Di^-^D'^dtik, (53) 

yh(L)-xhi^L) = 2D%, (54) 

yh{TV)i - xh(Tv)i = 6) (55) 

yh{TT)ij - xh(TT)ij = 0. (56) 

Since the gauge transformation rule (55) coincides with the gauge transformation rule (28) 
for the variable Xj, we may identify the variable Xi with h(Tv)i- 

Xi := h(Tv)i- (57) 

Thus, we have confirmed the existence of the variable Xi. Next, we show the existence of the 
variable Xt whose gauge-transformation rule is given by Eq. (27). Inspecting these gauge 
transformation rules (52) and (55), we find the definition of Xt as 

Xt := \vL)-^~^D^dt\ry)k- (58) 

Actually, the gauge transformation rule for Xt defined by Eq. (58) is given by Eq. (27). 
This is the desired property for the variable Xt. Thus, we have consistently confirmed the 
existence of the variables Xt and X^ which was assumed in the definitions (29)-(31) of the 
components of the tensor field Hah- This is the most non-trivial part of the ingredients of 
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this paper. Definitions (57) and (58) also imply that we may start the construction of the 
gauge-invariant and gauge-variant variables from the decompositions of the components hu 
and hij which are obtained by the substitution of Eqs. (58) and (57) into Eqs. (30) and (31) 
with the decomposition formulae (35) (37). This approach is discussed in See. 4. 

Now, we construct gauge-invariant variables for the linear-order metric perturbation hab- 
First, the gauge-transformation rule (56) shows that h(jT)ij is gauge invariant by itself and 
we define the gauge-invariant transverse-traceless tensor by 

Xij ■=h(TT)ij- (59) 

The transverse-traceless property of Xij is automatically given by the definition of h^x^^ij 
in Eqs. (36) and (37). Inspecting the gauge-transformation rules (53) and (55), we define a 
gauge-invariant vector mode Vi by 

J^i ■■= h(y)i-dth(Tv)i + DiA~^D''dth(Tv)k- (60) 

Actually, we can easily confirm that the variable is gauge invariant, i.e., yVi — xi'i = 0. 
Through the divergenceless property of the variable /i(y)i, we easily see the property DVj = 
0. The gauge-invariant variables for scalar modes are defined as follows: First, inspecting 
gauge-transformation rules (27) and (32), we define the scalar variable $ by 

-2$ := Htt - 2dtXt = Htt - '^dt\vL) + 2dtl^-^D^dt\TV)l- (61) 

Inspecting the gauge-transformation rules (54) and (55), we also define another gauge- 
invariant variable by 

-2n^ := h^L) - '^D'Xi = h^L) - '^D'h^TV)i- (62) 

We can easily confirm the gauge invariance of the variables $ and ^ through gauge- 
transformation rules (27), (32), (54), and (55). Here, we choose the factor of ^ in the 
definition (62) so that we may regard $ = ^' is Newton's gravitational potential in the 
four-dimensional Newtonian limit. 

In terms of these gauge-invariant variables and the variables Xt and Xj, which are defined 
by Eqs. (58) and (57), respectively, the original components {htt, hti, hij} of the metric 
perturbation hab is given by 

htt = -2^ + 2dtXt--{dta + pW,a- (3^/3'Ki,)Xt 

~ (^P'p'^P^Kkj - ma + aq'^dtPj 

+a^D'a - ap^D'Pk - P'P^DjOj X^, (63) 

hti = Ui + DiXt + dtXi--{Dia-P^Kij)Xt--M,^Xj, (64) 

a a 

hij = -2^qij + Xij + DiXj + D^X, + ^KijXt-^l3^KijXk. (65) 

Equations (63)-(65) imply that we may identify the components of the gauge-invariant 
variables T-Lab and the gauge-variant variable Xa so that 

Htt -.= -2^, nti-i^i, Hij := -2^qij + Xij (66) 
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and 

Xa:=Xt{dt)a + Xi{dx%. (67) 

These identifications lead to the decomposition formula (9) for the linear-order metric 
perturbation on an arbitrary background spacetime. 

We note that, in this outline of a proof, we assumed the existence of two Green function 
of the derivative operators A := D^Di and V^^ which is defined by Eq. (48). In other words, 
we have ignored the modes which belong to the kernel of these derivative operators A and 
V^K We call these modes as zero modes. To explicitly specify the Green functions for A and 
X'*-', we have to impose boundary conditions at boundaries (?E. Since the operators A and 
V"^ are elliptic, the change of the boundary conditions at dT, is adjusted by functions which 
belong to the kernel of the operators A and 2?*-^, i.e., zero modes. Therefore, we may say that 
the information for the boundary conditions for the Green functions A~^ and (P*-^) is also 
included in these zero modes. To take these modes into account, the different treatments 
will be necessary. We call this issue as zero-mode problem. 



3.2. Comparison with the FLRW background case 

Here, we consider the comparison with the ease where the background spacetime M.q is 
a homogeneous isotropic universe discussed in Refs. [9] to clarify differences of the above 
arguments in Sec. 3.1 from the well-known formulation of cosmological perturbations. The 
case of a homogeneous isotropic background universe corresponds to the case a = 1, = 0, 
and Kij = —Hqij, where H = dta/a and a is the scale factor of the universe. 

On this background spacetime, in Refs. [9], we decomposed the components hu and hij of 
the metric perturbation hab as 

hti = Dih{vL) + h(y)i, &h{v)i = 0, (68) 
hij = a'^h(^L)^ij + a'^h(T)ij, Y^h{T)ij = 0, (69) 

h(T)ij = (j^iDj - ^7ij^j \tl) + '^D^^jJurv)]) + ^(TT)ij, (70) 

D%TV)i = 0, D%TT)i3=^, (71) 

where Qij = a'^^ij, ^ij is the metric on a maximally symmetric space, Di is the covariant 
derivative associated with the metric 7jj, and A := D^Di. In the case where a = 1, = 0, 
and Kij = —Hqij, the decomposition (68)-(71) are equivalent to the decomposition (35)- 
(37) with Eqs. (29)-(31), (57), and (58) in this paper. Therefore, one might think that we 
may also apply the decomposition (68)-(71) even in the case of an arbitrary background 
spacetime. However, in the case of an arbitrary background spacetime, the decomposition 
(68)-(71) is ill-defined. Actually, if we regard that the decomposition (68)-(71) is that for an 
arbitrary background spacetime, we cannot separate /i(tl) and h{TV)j due to the non-trivial 
curvature terms of the background as pointed out by Dcscr [23]. These curvature terms 
come from the commutation relation between the covariant derivative Di and the derivative 
operator 2?*-^. This is why we apply the decomposition (35)-(37) with Eqs. (29)-(31), (57), 
and (58) instead of (68)-(71). 

Furthermore, in Refs. [9], we have assumed the existence of Green functions of the deriva- 
tive operators A, A -|- (n — 1)A;, and A -|- nA; to guarantee the one to one correspondence of 
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the set {hit, hu, hij} and {{hu, h^y^-j, h(^L), \tl)}, {h(y)u h(TV}i}, h(TT)ij}, where k is 
the curvature constant on the maximally symmetric space. The special modes which belong 
to the kernel of the derivative operators A, A. + {n — l)k, and A + nk were not included 
in the consideration in Refs. [9]. On the other hand, in this paper, we ignore the modes 
which belong to the kernel of the derivative operator A and P*-^. We note that the modes 
ignored in this paper coincides with the modes ignored in Ref. [9]. Trivially, the above oper- 
ator A := D^Di corresponds to the Laplacian A in this paper. In the case of the maximally 
symmetric n-space, the Riemann curvature and Ricci curvature are given by 

Rijki = 2kqi,[iqj]i = 2kqk[iqj]i, Rik = q^^Rijkl = {n - l)kqik. (72) 
In this case, the derivative operator V^^ defined by Eq. (48) is given by 

= (A + (n - l)k) +(^1-^^ D'D^. (73) 

When the operator V^^ acts on an arbitrary transverse vector field Vi {D^Vi = 0), we easily 
see that 

V'^Vj = ( A + (n - l)k) . (74) 

Therefore, the kernel of the derivative operator A + (n — l)k in Refs. [9] is included in the 
kernel of the derivative operator P*-' in this paper. Finally, we consider the case where the 
derivative operator V^^ acts on the gradient Dif of an arbitrary scalar function /: 



n 



D^A + ——R^^Di 

n — 1 



/• (75) 



In the case of maximally symmetric n-space, curvature tensors are given by Eqs. (72) and 
the derivative operator T>^^Di is given by 

V^^Dif = 2^^^^DUA + nk)f. (76) 
n 

To solve the equation V^^Dif = , we have to use the Green functions associated with 
the derivative operators A and A + nk. These are the reason for the fact that the Green 
functions A~^, (A -|- (n — 1)A;)~^, and (A -|- nk)^^ were necessary to guarantee the one-to- 
one correspondence between the components {hu, hij} and {/i(yz,), /j(y)i, /i(z,)) h^xL)^ h^xv)ij 
h(TT)ij} in Eqs. (68)-(71). In other words, we may say that the special modes belong to the 
kernel of the derivative operators A and V^^ are equivalent to the special modes which 
belong to the kernel of the derivative operators A, A -|- (n — l)k, and A -|- nA; in the case of 
the maximally symmetric space E in Refs. [9]. 

Finally, we note that the gauge-invariant variables defined by Eqs. (59)-(62) arc general- 
izations of the metric perturbation in the longitudinal gauge in cosmological perturbation 
theory [24]. As noted in Refs. [8, 16, 24], the choice of the gauge-invariant variable is not 
unique and there are many choice of the gauge-invariant variables. This situation corre- 
sponds to the fact that there are infinitely many complete gauge- fixing procedures. Actually, 
in Ref. [24] , through the specification of the gauge- variant parts AT" and Y"" of the first- and 
the second-order metric perturbations in Eqs. (9) and (15), we realized the two difi'erent 
complete gauge-fixing in the first- and the second-order perturbations, namely, the Poission 
gauge and the flat gauge, at kinematical level. Through the similar technique, Uggla and 
Wainwright [25] derived the linear-order Einstein equations both in the longitudinal gauge 
(Poisson gauge) and in the flat gauge (constant curvature gauge) in a compact manner. 
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4. Alternative construction of gauge-variant and gauge-invariant parts 

The result obtained in Sec. 3.1 implies that we may define the variables /i(yL), /j(y)i, ^(l)) 
^(Ty)i) ^(TT)ij by the following decomposition formulae for the components hti and hij: 



hti =: Dih(^vL) + h(^v) 
2 
a 



^ (Dia - p'^Kik) (/i(yL) " A-'D''dth^TV)k) " ^M, ''h^TV)k, (77) 



2 / \ 2 

+— -f^ij y^iVL) - At D dth(Tv)kJ - h{Tv)k, 

D^h(y)i = 0, q'^h(TT)ij = = D^h(TT)ij, 



(78) 
(79) 



where M--' is defined by Eq. (25). Here, these expressions are obtained through the substi- 
tution of Eqs. (57) and (58) into (30) and (31) and York's decomposition (35)-(37). In this 
section, we check the consistency of our result in Sec. 3.1. To do this, we change the starting 
point of our arguments to the decomposition formulae (77)-(79), though the starting point 
of our outline of a proof in Sec. 3.1 was Eqs. (29)-(31) with the assumption of the exis- 
tence of the variables and X^. First, we consider the derivation of the inverse relation of 
Eqs. (77)-(79) in Sec. 4.1. Then, we derive the gauge transformation rules for the variables 
h(VL), h(^Y)i, h(^L), h(TV)i, and h(TT)ij in Sec. 4.2. 

4.1. Inverse relation 

To derive the inverse relation of Eqs. (77)-(79), we first consider Eq. (77). Assuming the 
existence of the Green function for the elliptic derivative operator 



:= A - ^ (Aa - P^Kij) D' " 2£'' { ^ ( Aa - Kij) \ , 



(80) 



we obtain the relations 



1{VL) 



= 7-^ 



+A-'D%h^TV)k, 



hti - DiA-^D''dth^Tv)k + -Mi \TV)k 



+ 



Di-- [Dia - P^Kij) 



a 



a 



D^htk + D''dth^TV)k 



D^i-M^\TV)i 



a 



(81) 



(82) 



Equations (81) and (82) imply that we can obtain the relations between {/t(yL), and 
{hti, hij} if the relation between h^TV)i and {hti, hij} is specified. On the other hand, the 
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trace part and the traceless part of Eq. (78) are given by 



h{L) = q'^hij + -K/3''h(^TV}k 



--K ( F- 

a 



D'^htk - D%h(^rv)k + D'' ( -M^\tv)i 



1 2 ~ 



2 ~ 
a 



-1 



D'htk - D'^dth 



'tri(TV)k 



(83) 



(84) 



where we have used Eq. (81) and defined the traceless part Kij of the extrinsic curvature 
Kij by 



(85) 



Taking the divergence of Eq. (84), we obtain the single integro-differential equation for 
h{TV)k' 



V^%TV)k + D^ 



^Kmj I J'-'D'^ (^M^'h^TV)l - dth(TV)k 



-l3^h(Tv)k } 



\ fei 2 ~ 1 fe 

hmj (Irnjl h^l KjnjJ- D hfk 

n a 



(86) 



The existence and the uniqueness of the solution to this integro-differential equation 
is highly non-trivial. However, we assume the existence and the uniqueness of the solu- 
tion /icn/)^ = h(j'v)k [htm,hmn] to this integro-differential equation (86), here. This solution 

describes the expression of the variable h(Ty^i in terms of the original components hu and 
hij of the metric perturbation hab- Substituting the solution hfj'y^]^ = h^xy^}^ [htm,hmn] to 
Eq. (86) into Eqs. (81)-(83), we can obtain the representation of the variables h(yL-^, h(^v)ii 
/i(^) in terms of the original components hfi and hij of hab- Furthermore, the representa- 
tion of the variable /i(Tr)jj terms of hti and hij are derived from Eq. (84) through the 
substitution of the solution = ^{Ty)fc [htrn,h„in] to Eq. (86). 

Thus, the decomposition formulae (77)-(79) are invertiblc if the Green functions A~^, J-~^ 
exist and the solution to the integro-differential equation (86) exists and is unique. 



4-2. Gauge-transformation rules 

Through similar calculations to those in Sec. 4.1, we can derive the gauge-transformation 
rules for the variables /i(yL), h{L)i h{TV)ii h^^rp^ij. Prom Eqs. (81) and (82) the 

gauge-transformation rules (23) for the component hit, we obtain the gauge-transformation 
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rule for the variable h 



{VL) 



and h 



(V)i- 



yh{vL) 



yh(v)i 



a 



+ A-^D%Ak, 



(87) 



dtii - DiA-^D'^dtik 



+ 



A - - (Dia - P^Kij) 
a 



dtAk - -M^^Ai 
a 



-DiA-^D^dtAk 



a 



-M, ^Afc, 



where Aj 



yh(TV)i 



xh(TV)i ~ Ci- As in the case of the relations (81) and (82), these gauge- 
transformation rules (87) and (88) imply that we can obtain the gauge-transformation rules 
for the variable and /i(y)i if the gauge-transformation rule for the variable /i(Ty)i is 

specified. 

Prom Eq. (83) and the gauge-transformation rule (24), we can derive the gauge- 
transformation rule for the variable : 

2 



yhD-xh^L) = 2D'ii + -Kp^Ak 
2 



a 



+-K { J-^D'' 
a 



dtAk - -Mf}A, 



a 



(89) 



As in the case of the gauge-transformation rules (87) and (88), the gauge-transformation 
rule (89) also implies that we can obtain the gauge-transformation rule for the variable 

if the gauge-transformation rule for the variable /i(TV')j is specified. On the other hand, from 
the gauge-transformation rule for the traceless part (84) of hij, we obtain the equation 



+ yh{TT)ij - xh(TT)ij - -KijP^Ak 



2 ~ 



a 



dtAk - -Mk^Ak 



a 



0, 



where we have used Eqs. (23) and (24). The divergence of Eq. (90) yields 



a 



dtAk - -M^'Ak 

a 



0. 



(90) 



(91) 



Here, we note that we have assumed the existence and the uniqueness of the solution to 
Eq. (86). Since Eq. (91) is the homogeneous version of Eq. (86), this assumption yields that 
we have the unique solution = to Eq. (91), i.e., 

yh(TV)i - xh(Tv)i = ii- (92) 

Thus, we have specified the gauge-transformation rule for the variable h^xv)i ^^'^ the gauge- 
transformation rule (92) coincides with Eq. (49). 

Substituting Eq. (92) into Eqs. (87)— (90), we obtain the gauge-transformation rules for the 
variables /i(yi,), ^(L)) ^^'^ ^'(TT)ij- We easily see that the resulting gauge-transformation 

rules for these variables arc given by Eqs. (41), (42), (44), and (50), respectively. Further, 
we can construct gauge-variant and gauge-invariant parts of the metric perturbation hab in 
the same manner as in Sec. 3.1 and we can confirm Conjecture 1. Thus, we have reached to 
the same conclusions as those obtained in Sec. 3.1. Therefore, we may say that the results 
obtained in Sec. 3.1 are consistent. 
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5. Summary and discussions 

In summary, we proposed a scenario of a proof of Conjecture 1 for an arbitrary background 
spacetime which admits ADM decomposition. Conjecture 1 states that we already know the 
procedure to decompose the hnear-order metric perturbation hab into its gauge-invariant 
part T-Lab and gauge- variant part Xa- In the cosmological perturbation case, this conjecture 
was confirmed and then the second-order cosmological perturbation theory was developed in 
our series of papers [9, 16, 17]. However, as reviewed in Sec. 2, Conjecture 1 is the only non- 
trivial part when we consider the general framework of gauge-invariant perturbation theory 
on an arbitrary background spacetime. Although there will be many approaches to prove 
Conjecture 1, in this paper, we just proposed an outline a proof for an arbitrary background 
spacetime. 

In the outline shown in Sec. 3.1, we assumed the existence of Green functions of the elliptic 
derivative operators A and P*-'. This assumption implies that we have ignored the modes 
which belong to the kernel of these derivative operators. We call these modes as zero modes. 
To derive the explicit representation of these Green functions, we have to impose appropriate 
boundary conditions for the perturbative metric at the boundary OS. The modes which 
belong to the kernel of the above derivative operators also includes the degree of freedom 
of these boundary conditions, because of the cUipticity of the derivative operators A and 
V^^ as noted in Sec. 3.1. For this reason, we also emphasized the importance of these zero 
modes from a view point of the globalization of gauge- invariant variables [30]. Within the 
arguments in this paper, there is no information for the treatment of these mode. To discuss 
these modes, different treatments of perturbations will be necessary. We call this issue as 
zero-mode problem. The situation is similar to the cosmological perturbation case as noted 
in Sec. 3.2 and zero-mode problem exists even in the cosmological perturbation case. 

This zero-mode problem in cosmological perturbations also corresponds to the I = and 
1 = 1 (even) mode problem in perturbation theories on spherically symmetric background 
spacetimes. In the perturbation theory on spherically symmetric background spacetimes, 
one usually considers 2 -f 2 formulation [3] (or 2 -|- n formulation [4]), in which the similar 
decomposition to Eqs. (69)-(71) is applied. In this case, the indices ... in these equa- 
tions correspond to the indices of the components of a tensor field on 5^. Since S'^ is a 
2-dimensional maximally symmetric space with the positive curvature, we may regard n = 2 
and A; = 1 in Sec. 3.2. Then, the three derivative operators A, A -|- (n — l)k, and A -|- nA; in 
Sec. 3.2 are given by A, A -|- 1, and A -|- 2, respectively. Since the eigenvalue of the Laplacian 
A on 5^ is given hy A = + 1), we may say that the modes with I = and I = 1 belong 
to the kernel of the derivative operator A and A -|- nk, respectively. Therefore, we may say 
that the problem concerning about the modes with I = and / = 1 in the perturbations on 
spherically symmetric background spacetimes is the same problem as the zero-mode problem 
mentioned above. 

Thus, the arguments in this paper show that zero-mode problem generally appears in 
many perturbation theories in general relativity and we have seen that the appearance of 
this zero-mode problem from general point of view. Furthermore, as discussed in Appendix A, 
the kernel of the derivative operator T>^^ includes conformal Killing (or Killing) vectors on 
S. Therefore, the existence of zero mode is also related to the symmetry of the background 
spacetime. To resolve this zero-mode problem, careful discussions on domains of functions for 
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perturbations and its boundary conditions at dT, will be necessary. We leave this zero-mode 
problem as a future work. 

In our outline of a proof in Sec. 3.1, we used a tricky logic to construct gauge-invariant 
and gauge- variant parts of the metric perturbation hab- Therefore, in Sec. 4, we confirmed 
the consistency of our result in Sec. 3.1 through the change of the starting point of our 
arguments, and then, we have reached to the same conclusion as in Sec. 3.1. Due to this 
fact, we may say that the result obtained in Sec. 3.1 is consistent. In the approach in Sec. 4, we 
assume the existence of Green functions A""^ and T"^. Further, we also assume the existence 
and the uniqueness of the solution to the integro-differential equation (86). Although the 
correspondence between the existence of the Green functions which are assumed in Sec. 3.1 
and these assumptions is not clear, we may say that the above zero-mode problem is delicate 
in the case of an arbitrary background spacetime. 

We also emphasize that the gauge issue discussed in this paper is just a kinematical one 
and we did not use any information of the field equations such as Einstein equations. In 
other words, the results in this paper is also applicable not only to general relativity but 
also to any other metric theories of gravity with general covariance. 

Since our motivation of this paper is in higher-order general-relativistic perturbation the- 
ory, readers might think that the ingredients of this paper is related to the issue so-called 
"linearization instability" [27]. However, at least within the level of the ingredients of this 
paper, our arguments themselves still have nothing to do with the issue of linearization 
instability. The linearization instability is the issue of the existence of the solutions to the 
initial value constraints in the Einstein equations, while our arguments in this paper is just a 
kinematical one and have nothing to do with field equations as mentioned above. Although 
it will be interesting to reconsider the issue of linearization instability from the view point of 
our gauge-invariant perturbation theory, this issue is beyond the current scope of this paper. 
To discuss the linearization instability from our view point, it will be necessary to discuss 
the above zero-mode problem, at first. 

Although we should take care of the zero-mode problem, we have almost completed the 
general framework of the general-relativistic higher-order gauge-invariant perturbation the- 
ory. The outline of a proof of Conjecture 1 shown in this paper gives rise to the possibility of 
the application of our general framework for the higher-order gauge-invariant perturbation 
theory not only to cosmological perturbations [9, 16, 17] but also to perturbations of black 
hole spacetimes or perturbations of general relativistic stars through an unified formulation 
of the general-relativistic perturbation theory. Furthermore, as mentioned above, the results 
in this paper is also applicable not only to general relativity but also to any other metric 
theories of gravity with general covariance. Therefore, we may say that the wide applica- 
tions of our gauge-invariant perturbation theory will be opened due to the discussions in this 
paper. We also leave these development of gauge-invariant perturbation theories for these 
background spacetimes as future works. 
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A. Covariant orthogonal decomposition of symmetric tensors 

Since the each order metric perturbation is regarded as a symmetric tensor on the background 
spacetime {^4o,gab) through an appropriate gauge choice, the covariant decomposition of 
symmetric tensors is useful and actually used in the main text. Here, we review of the 
covariant decomposition of symmetric tensors of the second rank on an curved Riemannian 
manifold based on the work by York [26]. 

On an arbitrary curved Riemannian space {T,,qab) (dimS = n), one can decompose an 
arbitrary vector or one-form into its transverse and longitudinal parts as 

Aa = Ai{dx')a = {D.A^L) + ^(y)i) {dx\, D'A^v)i = 0, (Al) 

where Di is the covariant derivative associated with the metric qab = qij{dx^)a{dx^)h- ^(l) 
is called the longitudinal part or the scalar part and Afy^i is called the transverse part or 
vector part of the vector field Aa on (Ti,qab)-, respectively. 

Moreover, this decomposition is not only covariant with respect to arbitrary coordinate 
transformations, it is also orthogonal in the natural global scalar product. To clarify this 
orthogonality, York [26] introduced the inner product for the vector fields on S. This is, for 
any two vectors V"- and W"-, we have 

{V,W) := J^e^V^W^b, (A2) 

where eg denotes the volume element which makes the integral invariant and the integration 
extends over the entire manifold (S,^^^). In terms of this inner product, the orthogonality 
of the vector fields V = D^V^l) ■= <f^DbV^L) and W = q"-^V(y)b with D''V(y)a = is given 
by 

/ €qDaV^L)V(y)bq''^ = [ SaV^L)V^v)bq'"' - [ e^V^L)DaViy)m''\ (A3) 

where is the volume element of the (n — l)-dimensional boundary OS of S. Since the 
second term of Eq. (A3) vanishes due to the condition D'^V^y^a = 0) the inner product (V, W) 
vanishes if V(2.) and V(y)b satisfy some appropriate boundary conditions at the boundary 9S 
of S so that the first term of Eq. (A3) vanishes. In this sense, the scalar part (the first 
term in Eq. (Al)) and the vector part (the second term in Eq. (Al)) orthogonal to each 
other. Geometrically, the decomposition of 1-forms, and more generally p-forms, leads via dc 
Rham's theorem to a characterization of topological invariants of S (i.e., Betti Numbers) [28]. 

In this appendix, it is assumed that the n-dimensional space S is closed (compact manifolds 
without boundary) following York's discussions. Although the decomposition discussed here 
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will also be valid for the other n-dimensional spaces S with the boundary dT, with some 
appropriate boundary conditions at dH, in this Appendix, we choose the closed spaces as the 
topology of S for mathematical convenience. Through this assumption, in this Appendix, we 
consider the TT-decomposition (transverse traceless decomposition) of a symmetric tensor 
tp"''^ on S, which is defined by 

V'"' = <T + < + (A4) 

where the longitudinal part is 

n 

and the trace part is 

^ab^.= l^^qab^ ^..= q^^^cd. (A6) 

Let us suppose that both an arbitrary symmetric tensor field V'"'' and the metric Qab are 
C°° tensor fields on S. First, we define ip^j, in accordance with Eq. (A4) by 

^l,f^ ■= - 1^5"^ - {LWfK (A7) 
We note that the tensor ij^^j, is traceless, i.e., 

Qabi^fr = (A8) 
by its construction (A7). Further, we require the transversality on the tensor field V'tT' ^■^■■> 

DbtpfT = 0- (A9) 
Equation (A9) leads to a covariant equation of the vector field VF" in Eq. (A7) as 

DaiLWr' = Da - ^V'?"') • (AlO) 



The explicit expression of (AlO) is given by 

V'^'^Wc = Da i^r'' - , (All) 

where the derivative operator V^'^ is defined by 

where i?^^ is the Ricci curvature on (E, §„;,). 

The basic properties of Eq. (All) arc also discussed by York [26]. The operator 2?°^ defined 
by Eq. (A12) is linear and second order by its definition. Further, this operator is strongly 
elliptic, negative-definite, self-adjoint, and its "harmonic" functions are always orthogonal to 
the source (right-hand side) in Eq. (All). Here, "harmonic" functions of V"-^ means functions 
which belong the kernel of the operator V"-''. Moreover, he showed that Eq. (All) will always 
possess solutions W"' which is unique up to conformal Killing vectors. Due to this situation, 
in this paper, we assume that the Green function {'D~^)ab defined by 

{V-\,V^ = V'\V-\, = 5a'' (A13) 

exists through appropriate boundary conditions at the boundary 9S of S. Although York's 
discussions are for the case of the closed space S, we review his discussions here. In this 
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review, we explicitly write the boundary terms which are neglected by the closed boundary 
condition to keep the extendibility to non-closed S case of discussions in our mind. 

The ellipticity of an operator depends only upon its principal part, i.e., the highest deriva- 
tives acting on the unknown quantities which it contains. To see the ellipticity of an operator, 
we consider the replacement of the each derivative operator Da occurring in its princi- 
pal part by an arbitrary vector V^. Through this replacement, the principal part of the 
operator defines a linear transformation cr„. The operator is said to be elliptic if (T„ is an 
isomorphism [29]. In the present case, 

[a^.iV)]"'' = V^V'' + q^^VcV''. (A14) 

Here, operates on vector Xa and defines a vector-space isomorphism when the determinant 
of (7„ is non- vanishing for all non- vanishing V"'. We can verify deta^ / 0, for example, by 
choosing = {d/dx'^)"' in a local Cartesian frame {x'^}. The operator is said to be strongly 
elliptic if all the eigenvalues of are nonvanishing and have the same sign. This is easily 
checked and 'D"'^ is strongly elliptic. 

To show that D"'^ is negative definite, we consider the inner product (A2) of the vector 
field VW" := V^Wb and 

= [ ScWb{LWf'-l [ egiLW)ULWf', (A15) 

where we use the fact that the tensor (LW)'"^ is symmetric and traceless. Eq. (A15) shows 
that the operator V"'^ has the negative eigenvalues in the case where the first term (boundary 
term) in Eq. (A15) is neglected, unless (LW)^'^ = 0. The self-adjointness of the operator V"'^ 
is follows from a similar argument in which one integrates by parts twice: 

f e,qabV%VW)' = [ SchiLW)''^-iLV)''Wb 

+ J^egW^V'^Vc (A16) 

for any vectors V and W, where we use the fact that the tensor (LW)"'^ and [LV)"'^ arc 
symmetric and traceless. Eq. (A16) shows that the operator is self-adjoint if the first 
term (boundary term) in Eq. (A16) is negligible. 

When we can neglect the boundary terms in Eq. (A15), the right-hand side of (A15) can 
vanish only if {LW)"'' = 0. This means either W"" = or is a conformal Killing vector (or 
Killing vector) of the metric qab- The condition for a conformal Killing vector is, of course, 
not satisfied for an arbitrary metric but this is given by 

£wqab = Mab (A17) 

for some scalar function A, where £w denotes the Lie derivative with respect to W"-. Taking 
the trace of both sides, we find 

A = -DcW. (A18) 

n 

Therefore, W"' is a conformal Killing vector if and only if 

D^'W'' + D^W - '^q^-^DcW = {LWy^ = 0. (A19) 
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It follows that the only nontrivial solutions of V^^Wb = are conformal Killing vectors if 
they exist. Hence the nontrivial "harmonic" functions of T'"'^ are conformal Killing vectors. 

We shall now show that even if these "harmonic" solutions exist, they are always orthogonal 
to the right-hand side of (AlO) and, hence, can cause no difficulties in solving equation (AlO) 

by an eigen function expansion. Denote the conformal Killing vectors by W"" = C", where 
by definition {LCY^ = 0. Form the scalar product of the right-hand side of (AlO) with C 
and integrate by parts to find 

= 0, (A20) 

where we use the fact that V'"* ~" ^Q^'^fp is symmetric and traceless and we also neglect the 
boundary term. Hence the source in Eq. (All) is in the domain of {T>~^)"''' and (X>^^)"^ gives 
the solution to Eq. (All) even in the presence of conformal symmetries. 

These results also show that the solution to Eq. (All) must be unique up to conformal 
Killing vector fields. Since only (LW)"''^ enters in the definition (A7) of V'tT' conformal 
Killing vectors cannot affect ip^j,. 

The orthogonality of V'tt' (LW)'^'', and ^V'^"^ is easily demonstrated. We see readily that 
^tpq""^ is pointwise orthogonal to (LW)"''' and to ^'tt' ^ (LW)^^ and V'tt both trace- 
free. To show that i/^jij, and (LF)"'' are orthogonal for any vector V"' and any TT tensor, 
we only show that 

= 0, (A21) 

where we use the fact that the tensor -0^^ is symmetric, traceless, and transverse (A9). We 
also neglect the boundary term in Eq. (A21). Thus, we conclude that the decomposition 
defined by (A7) exists, is unique, and is orthogonal. 

One can further decompose the vector W"" uniquely into its transverse and longitudinal 
parts with respect to the metric Qab- This splitting is orthogonal, as in Eq. (Al). 

Since the above discussions are for the case of a closed spaces S, careful discussions on the 
boundary terms, which are neglected in the case of a closed S, is necessary if we extend the 
above arguments to the case of a non-closed S. However, we do not go into these detailed 
issues. Instead, in the main text, we assume that the existence of the Green function of 
the derivative operator and use the transverse-traceless decomposition for an arbitrary 
symmetric tensor on S discussed here. 
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